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fields,  the  ealibniion  curves  and  derivative  methods  for  the  caae  ^  in  magnetic 

The  advantages  and  typical  error,  of  both  technique  5  nn^g  m^etograph. 
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of  the  SVMG.  ^  present  instrumental  configuration 


INTRODUCTION 

Vector  magnetic  fields  are  a  leitmotiv  in  modem  solar  physics  esDcciallv  for  their  ml.  •  i 
olution  (Zim,  1985;  Ai,  1987  Hxgyard,  igsTRiefS 

However,  two  issues  are  still  the  nightmare  of  every  physicist  that  wantu  to  f?..i  -ti.  * 
i,  k.  (.k.)  i.  I  ‘ 

1-  the  n^  to  get  observations  of  extremely  high  quality  an^ 

Sxc  PUk.  U.,  c<,n.id.ri.g  iU  •' 

works  at.  Nevertheless  we  hone  that  anm.  of  th  -j  ?  spectral  line  the  instrument 

«.efk.  fo,  oth,,  *"-•  k=.  mad.  .ill  ,„;d, 

.tc.)  kx^xlrepd,  “‘’.rSi*'  “f  k™ 

mque.,  Jnl  for  ”lk«>reti»l"  condiliom,  tk.i.  cojSdeiiii.  ihl  ^Ztlr^ 

error  someth.  ®  *  effect  of  spectral  smearing  and  other 
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1.  THE  INSTRUMENT  AND  THE  CALIBRATION  METHODS 


The  JHU-APL  Sol*r  Vector  M&gnetograph  (SVMG)  U  an  imaging  magnetograph,  i.e.  it  ac- 
quires  nanow  band  images  in  different  Ught  polarisation  states  (selected  through  different 
of  a  quarter  waveplate  and  a  Glan-Taylor  prism).  The  images  correspond  to  measmementi  of  I+Q, 
I-Q  I+U  I-U,  I+V,  I-V,  where  I,  Q,  U,  V  are  the  usual  Stokes  parameters,  and  are  acquired  at 
a  rate  of"'-  30  images/min.  In  a  normal  observing  procedure,  up  to  20  images  are  accumulated 
for  each  polarisation  state,  in  order  to  increase  the  signal-tc^noise  ratio;  this  brings  the  acquisition 
time  for  a  sequence  to  ~  5  minutes.  The  field  of  view  (FOV)  u  ~  4  x  3  aremin’.,  the  pixel  sue  is 


120  mA 


0.4arc8ec.  .  .  t 

The  SVMG  operates  in  the  magnetic  Cal  6122.22  A  line,  that  has  a  Doppler  width  of  ‘ 

and  a  Land^  factor  jx,  of  1.75.  ^  ' 

The  spectral  selection  in  the  SVMG  U  performed  by  a  tunable,  bthium-mobate,  Fabry-Perot  In¬ 
terferometer  (FPI)  and  by  an  interference  blocker,  that  acts  as  a  prefilter  for  the  FPL 
The  FPI  has  a  transmission  profile  described  by  the  usual  Airy  function  (Hernandei,  1986),  with  a 
FWHM  of  ~  150  mA  and  a  FSR  of  3.25  A  (Fig.  la).  The  transmission  is  nominally  around  50%, 
but  no  precise  measurements  have  been  recently  made. 

Duncan  and  Thomas,  1987,  provide  some  information  about  the  spectral  position  of  the  FPI  orders 
in  the  range  6000-7000  A. 


The  blocker  has  a  FWHM  of  ~  lA  and  a  FSR  of  16  A  (Fig.  lb).  The  only  ways  to  tune  the 

blocker  are:  .  ,  i /j  \ 

1-  tilt  it,  with  the  effects  of  shifting  the  passband  towards  the  blue  (~  20  mA/deg),  or 

2-  changing  the  temperature,  moving  the  passbands  of  approximately  80  mA/*C  (cooling  -  blue 


shift;  heating  -  red  shift),  ,  ,  i,  .  ..  u 

The  second  option  is  not  feasible  during  an  observing  sequence,  and  the  first  hasn  t  been  mves^ 
gated  carefully  up  to  now,  so  that  actually  during  the  observations  the  blocker  sits  at  one  fixed 
wavelength.  The  central  position  of  the  blocker  passband  is,  as  we  will  see  later,  important  for  the 
validity  of  the  calibration  procedure,  and  it  must  be  carefully  chosen. 


The  nature  of  the  SVMG,  an  imaging  filter  magnetograph,  doesn’t  allow  the  usei^f  calibration 
methods  based  on  a  full  profile  analysis  of  the  Stokes  profiles  I,  Q,  D,  V  (see  Skumanich  and  Lites, 
1987),  at  least  on  a  regular  basis,  for  two  main  reasons: 

1-  the  acquisition  of  images  along  the  line  profile  with  a  good  spectral  sampling  rate  1/2  FWHM, 
50-70  mA)  will  take  too  much  time  (more  than  1  hr.),  with  all  the  problems  related  to  changes  of 


the  structures 

2-  the  computing  time  necessary  to  perform 
the  FOV  would  be  a  prohibitive  restriction. 


this  analysis  for  all  the  2  x  10*  (or  more)  points  in 


The  daily  calibration  must  be,  therefore,  based  on  the  signal  at  1  or  few  wavelengths  m  particular 
positions  in  the  line. 

Two  possible  methods  are: 

1-  conversion  of  polarisation  signal  to  magnetic  field  through  calihntion  curves,  built  solving 
the  equations  of  transfer  for  polarised  radiation  for  the  spectral  line  in  a  particular  atmospheric 
model,  and  finally  considering  the  spectral  transmission  profile  of  the  instrument  (in  our  case  FPI 
+  blocker).  This  is  the  method  adopted,  for  example,  by  the  Marshall  group  (Hagyard  et  al., 
1988),  and,  for  the  limited  case  of  longitudinal  field,  by  many  others  (Title  ei  al.,  1987;  Lundstedt 

ei  al.,  1991). 

2-  in  weak  field  approximation  (WFA),  use  of  the  derivative  method,  that  states  the  proportionality 
of  the  polarisation  signals,  both  linear  and  circular,  to  magnetic  field,  through  the  slope  of  the  line 


‘An  slteniaUve  "iolutloo”  could  be  the  acquiution  of  imaje*  alon*  tbe  whole  line  profile  with  a  imaller  nui^er 
of  accimmlationj  (even  tinaJe  image*),  that  will  reduce  the  total  time  to  ~  10  min.  According  to  Ute*  ^d^- 
manich  (1984),  thi*  will  »till  give  a  precition  of  5%  (for  field*  bigger  than  1000  G)  and  if*  »omething  worth  further 
inve*tigsUon*. 
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profile  (JefTenea  et  •/.,  1989;  Jefferies  &nd  Mickey,  1991,  bere&fler  JM). 

UnfortuBately,  at  this  time  the  SVMG  is  still  ib  the  "test*  phase  so  that  few  observations  are  bow 
available.  In  particular,  no  observations  along  the  whole  line  profile  have  been  performed,  and  a 
check  of  the  validity  of  the  2  mentioned  calibration  methods  vs.  fitting  line  profiles  methods  is  not 
possible  yet. 


2.  SYNTHETIC  LINE  PROFILES 

The  profiles  have  been  computed  using  Graham  Murphy’s  code,  Stokes  Profiles  Synthesis  Rou¬ 
tine  (SPSR,  Rees  et  of.,  1989),  for  the  Cal  6122  line  and  three  different  atmospheric  models:  quiet 
sun  VAL3C  (Vernassa  et  of.,  1981),  Moe-Maltby  penumbral  (1969)  and  Haltby  umbral  (Avrett  et 
aJ.,  1986).' The  magnetic  parameters  have  been  assumed  constant  with  height  in  all  the  calcula¬ 
tions;  we  used  a  grid  of  values  of  B  (total  magnetic  field)  and  j  (mclination  of  magnetic  field  with 
respect  to  the  line  of  sight)  ranging  from  0  to  2500  G  (step  50  G),  and  from  0  to  90*  (step  30®) 
respectively.  The  asimuth  of  the  field,  x,  i«  normally  computed  through  the  ratio  of  the  parameters 
Q  and  U  independently  from  the  adopted  calibration  technique,  and  therefore  we  limited  ourselves 
to  the  case  of  x  =  0®.  Nevertheless,  we’ll  make  a  brief  comment  about  the  determination  of  this 
parameter  and  errors  involved,  in  Sect.  5. 

The  wavelengths  step  for  the  simulated  profiles  is  10  mA. 

Fig.  2a  shows  the  I  profile  for  the  VAL3C  model,  and  no  magnetic  fields;  Pigs.  2b-d  show  Stokes 
I  and  V  (on  the  same  scale)  for  3  different  seta  of  B,  y  and  x-  Figs.  3b-d  show  the  same,  for  the 
U  and  Q  parameters.  All  the  profiles  are  normalised  to  the  contmuum  intensity. 

It  might  be  useful  to  remind  that  for  7  =  0®  the  transverse  field  is  absent,  while  when  y  =  90®  the 
longitudbal  field  is  sero. 

Figs.  4-5  and  6-7  are  the  same  as  2  and  3,  but  for  the  Moe-Maltby  penumbral  and  Maltby 
umbral  models,  respectively.  The  line  is  enormously  strengthened  in  the  penumbral  and  umbral 
models,  respect  to  the  quiet  sun,  and  this  can  be  a  font  of  errors  in  the  calibration,  especially  using 
the  calibration  curve  method.  The  validity  of  umbral  models  is  somehow  in  discussion,  but  now 
we  don’t  have  any  observation  of  the  line  profile  in  a  spot,  that  can  confirm  or  discharge  these 
simulations. 

3.  CALIBRATION  CURVES.  POLARIZATION  VS.  MAGNETIC  FIELD. 

3.a  The  ’’theory” 

The  calibration  curvet  method  essentiaUy  consists  in  determining  functional  relationships  be¬ 
tween  the  magnetic  field  and  the  correspondent  fractional  polarisation  (both  circular  and  linear). 
The^  relationships  are  derived  solving  the  equations  of  transfer  for  polarised  radiation,  for  the 
particular  case  of  Constance  of  magnetic  parameters  with  height  in  the  atmosphere,  as  mentioned 
in  Sect.  2.  We  can  write: 

Bl  =  fi{y.  A,  B,  ..)Py  (1) 

Bt  =  My,  A,  B,  ..)P^/*  (2) 

where  JV  =  V/I,  is  the  fractional  circular  polarisation,  Pq  =  y/Q^  +  U^/I  is  the  fractional  linear 
polarisation,  and  Bi  and  Br  the  longitudinal  and  transverse  magnetic  fields. 

/i  and  A  are,  in  general,  functions  of  the  atmospheric  model,  the  magnetic  parameters  (for  example, 
iaere  is  a  st  ng  dependence  on  the  mclination  7),  the  transmission  profile  of  the  filter  and  the 
wavelength  position  m  the  spectral  line.  Examples  of  such  curves  are  shown  in  Figs.  8  and  9,  for 
the  Moe-Maltby  penumbral  model,  and  Bj,  and  Br  respectively.  The  fractional  polarisation  is 
shown  for  several  spectral  positions  in  the  blue  wing  of  the  line  (from  -260  mA  to  -120  mA  from 
line  center);  the  spectral  smearing  due  to  the  FPI  and  blocker  hasn’t  been  taken  into  account. 
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Pig.  10  is  the  same  as  Fig.  0,  but  on  a  logarithnric  scale,  to  better  show  the  (quite  smaU)  values 
of  linear  polarisation. 

In  weak  field  approximation  (WFA)  Eqs.  (1)  and  (2)  become  simpler  (Hagyard  tt  ai.,  1988): 

Bi  =  d/V  (3) 

Bt=C2P;/^  (4) 

where  Ci  and  cj  are  now  constant  with  the  magnetic  parameters,  even  though  still  dependent  on 
the  wavelength  and  the  filter  characteristics  (and  the  model,  as  before).  The  WFA  nominally  holds 
for  values  of  vb  —  '<^1|  where  is  the  splitting  Zeeman,  or 

AAb  =  mA»px,B  =  4.67  X  10-‘»AVB  ,  (5) 

with  A  measured  m  A  and  B  in  Gauss’.  For  the  Ca  6122  line,  this  translates  b  B<  2000  G, 
assu^g  a  value  of  120  mA  for  AAp.  However,  a  quick  look  at  Figs.  8  and  10  will  tell  us  that 
this  is  not  strictly  true,  i.e.  the  linearity  between  Bi  and  the  circular  polarisation  holds  up  to  - 
100(^1500  G,  but,  even  b  this  range,  the  differences  due  to  the  bclbation  are  relevant.  The  same 
considerations  apply  to  the  ,  where  the  deviations  from  the  Ibearity  are  more  severe. 

Usually  the  curves  of  polarisation  vs.  magnetic  field  are  shown  for  just  one  value  of  bclination 
(7  =  0®  for  the  longitudinal  field,  and  7  =  90®  for  the  transverse),  and  no  mention  is  made  about 
the  eventual  differences,  but  this  is  a  source  of  significant  error,  as  we  will  see  b  Sect.  3.d. 

Some  prelimbary  tests  showed  that  the  same  problems  are  present  b  the  magnetic  Ibe  Pel  5250 
also,  so  we  can  probably  generalise  the  question  of  the  validity  of  relations  (3)  and  (4). 

The  choice  of  the  working  wavelength  is  dictated  by  the  need  of  a  good  signal- to-noise  ratio,  and  by 
the  need  of  avoiding,  as  much  as  possible,  saturation  effects,  magneto-optical  effects,  dependence 
on  bclination  values  etc;  its  value  must  be  carefully  considered  b  order  to  satisfy  these  conditions. 
In  Figs.  8  and  9  is  clear  that  the  signal  is  bereasbg  noticeably  gobg  towards  Ibe  center,  but  at 
the  same  time  some  of  the  problems  just  mentioned  appear:  the  difference  b  the  curves  due  to 
the  7  values,  for  example,  is  bereasbg  with  wavelength,  especially  for  the  longitijdbal  field;  the 
magneto-optical  effects  are  present  at  A=6122.08  and  6122.1  A,  etc.  ^ 

Due  to  the  particular  combination  of  a  big  AA^j  and  a  quite  low  gi  factor,  th^Tcal  6122  Ibe 
doesn’t  suffer  of  saturation  problems,  unless  we  go  too  close  to  the  Ibe  center  br  at  very  high 
magnetic  fields  (>  3000  G).  The  same  two  factors,  however,  contribute  to  rather  small  values  of 
polarisation:  b  Fig.  8  we  can  sec  that,  at  A=6122.02  A,  for  example,  we  have  a  V/I  ratio  of  ~ 
3%  for  500  Gauss;  for  the  magnetic  Fcl  5250  line  (yi=3  and  AA^,  -  40  mA),  at  a  similar  distance 
from  the  Ibe  center,  the  fractional  circular  polarisation  would  be  close  to  15%  for  the  same  500  G. 
The  values  of  Ibear  polarisation  are  better  visualised  b  Fig.  10,  and  are  not  very  encouragbg:  at 
A=:6122.02  A,  a  sensitivity  of  lO’  will  be  sufficient  to  detect  transverse  fields  of  no  less  than  3-400 
G. 

The  FPI  and  the  blocker,  though,  change  the  Stokes  parameters  profiles:  we  must  hence  analyse 
their  effects,  and  choose  our  workbg  wavelcngth(B)  on  the  base  of  the  results. 

We  will  analyse  the  effects  of  the  FPI  alone,  first,  and  then  add  the  blocker.  This  could  be  useful 
b  a  near  future,  when  a  tunable  blocker  will  be  allowable,  and  the  net  spectral  transmission  profile 
will  be  the  one  of  the  FPI  only. 

3.b  Effects  of  spectral  smearing.  FPI 

To  simulate  the  effects  of  the  spectral  smearbg  btroduced  by  the  Fabry-Perot,  we  convolved  the 
Stokes  parameters  correspondent  to  the  different  models  and  values  of  magnetic  fields,  described 
b  Sect.  2,  with  the  transmission  profile  shown  b  Fig.  la.  It  might  be  useful  to  rembd  that  we 
must  convolve  the  parameters  and  then  do  the  ratio  to  obtab  the  fractional  polarisations  of 

*FoDoiraj  Jefferiei  tt  *1,  1989,  0.5  is  a  reasonable  upper  limit  for  ifg. 
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Figt.  8-10:  thU  is  because  the  ratio  and  the  convolution  do  not  commute. 

Fig.  11  shows  the  Stokes  profiles  for  the  case  of  B=1500  G,  7  =  30*  and  x  =  W*,  in  photospheric 
model,  "original*  and  convolved  with  an  Airy  function  of  150  mA  of  FWHM  (and  normalising 
the  continuum  intensity  to  1,  as  before).  The  profiles  are  shallower  as  expected,  but,  happily,  the 
convolved  V,  U,  Q  signals  are  greater  than  the  original  ones  for  wavelengths  <  8122.06  A  (this 
•cutting"  wavelength  depends  on  the  model  and  the  filter  characteristics);  in  threase  of  U  and  Q, 
the  smearing  profile  is  big  enough  to  dilute  all  the  details  m  the  line  cent^. 

Figs.  12-14  are  the  analogous  of  Figs.  8-10,  this  time  considering  the  convolution  with  the  PPI 
profile.  The  polarisation  signal  is  in  this  case  bigger  for  wavelengths  <  6122.04  A.  We  can  notice 
from  Pig.  12  that  the  range  of  linearity  for  is  greatly  extended  (up  to  the  limit  of  our  simulations, 
i.e.  2500  G),  and  the  differences  due  to  the  7  vsJues  are  reduced,  even  if  not  completely.  Por  the 
transverse  field,  the  same  extension  of  the  linear  regime  applies,  but  the  differences  introduced  by 
various  7  values  are  enhanced;  the  reason  for  this  apparently  strange  behavior  sits  in  the  rapidly 
varying  shape  of  the  U  and  Q  profiles,  and  in  their  dependence  on  the  bclination. 

We  can  probably  clarify  the  issue  with  an  example:  Pig.  15a  shows  two  Q  profiles  for  Br=1000 
G,  7  =  30*  and  90®,  and  X  =  0*-  The  two  profiles  are  extremely  different  but  m  the  line  center 
and  the  far  wings,  where  we  ought  work.  However,  when  we  convolve  with  the  PPI  profile,  whose 
FWHM  is  bigger  or  comparable  to  the  wavelength  scale  of  variation  of  the  Q  parameters,  the  two 
signal  do  not  coincide  anymore,  even  in  the  far  wings,  as  shown  in  Pig.  15b^.  Since  we  don’t 
have  any  other  independent  way  to  determine  7,  we  must  choose  one  or  another  of  these  curves 
(or  a  combination  of  them)  as  our  ealihniion  curve:  this  will  lead  us  in  any  case  to  errors  in  the 
determination  of  magnetic  fields. 

The  spectral  smearing  hence,  as  we  will  see  in  Sect.  4  also,  has  double-faced  effects:  on  the  good 
side,  it  enhances  the  signal  level  in  the  wings  and  increases  the  linearity  of  the  curves,  therefore 
reducing  errors  due  to  the  velocity,  e.g.;  on  the  bad  side,  it  brings  signal  from  the  line  center, 
with  all  the  related  problems,  to  the  line  wings,  therefore  introducing  dishomogeneities  like  the  7 
dependence. 

3.C  Effects  of  spectral  smearing.  Blocker 

In  the  actual  version  of  the  SVMG  the  blocker  is  not  easily  tunable,  so  we  have  to  determine 
which  is  the  most  convenient  wavelength  to  sit  the  blocker  on.  Fig.  16  shows  the  same  case  of  Fig. 
11,  for  two  different  positions  of  the  blocker,  m  line  center  and  in  the  blue  wing  (note  the  different 
scale  on  the  y  axis  respect  to  Fig.  11).  For  the  intensity  profile  the  differences  are  15-30%  when  the 
blocker  is  position  at  line  center  (in  the  range  of  wavelengths  6121.96-6122.10),  and  5-10%  when  m 
the  wing;  the  differences  in  V,  Q,  U  are  not  as  big  (5-10%  of  the  original  polarisation  signal  when 
the  blocker  is  positioned  at  line  center,  and  2-5%  when  in  the  wing).  However,  the  quantities  Py 
or  Pq  will  be  different  from  the  case  of  FPI  only,  so  the  calibration  curves  must  be  built 
into  account  the  blocker. 

If  we  work  at  one  single  wavelength,  or  few  of  them  very  close  together,  the  best  position  for  the 
blocker  is  that  same  working  wavelength;  if  we  want  to  perform  a  scan  of  the  whole  line  profile,  it 
will  be  easier  to  keep  the  blocker  at  line  center,  in  order  to  maintain  things  symmetric. 

Figs.  17-19  are  the  same  as  Figs.  12-14,  but  taking  into  account  the  influence  of  the  blocker,  set 
at  the  same  nominal  wavelengths;  apart  for  slight  differences  in  the  fractional  polarisation  signals, 
the  curves  have  the  same  shapes. 


In  an  the  umulsUonj,  the  contrwt  of  the  intenuty  profile  (defined  ai  -  /,  *  )//««».)  i*  around  35%, 

whAe  the  ob.er»ed  one  doean’t  go  over  a  mere  15%.  Thii  luggeati  the  preaence  of  a  itrong  percentage  of  icattcted 
light,  that  will  affect  all  the  calibration  procedures.  Ifa  therefore  crucial  try  to  reduce  this  scattered  light  at  the 
minimum  possible  level  by  means,  e.g.,  of  several  field  stops  in  the  optical  path.) 

‘This  problem  is  common  to  all  the  imaging  magnetographs,  where  the  spectral  selection  is  performed  hr  means 
of  filters  with  pastbands  around  100  mA. 
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S.d  Calibration  curves  and  errors 

.***  I*"'  wavelengths  are  m  the  range  6122.00- 

61M.06  A,  in  bet,  at  lower  wavelengths  the  signal  u  extremely  smaU,  especially  for  the  linear 

^  j  ^  -350  G.),  while  cl^  to  the  line  center 

the  inclination,  magneto-optical  and  saturation  effects  are  stronger 

Sfl  iMlruma.t  «  ~  2  -  3  x  l(l>,  .t  tka.  .h„„ld  fc. 

able  to  detect  transverse  fields  of  1-200  G,  and  longitudinal  fields  of -n.  lo  G. 

^bl«  1  and  2  show  the  values  of  the  calibration  constants  cj  and  cj  of  Eqs.  (3)  and  (4)  ob¬ 
tain^  from  our  simulations.  Some  more  or  less  obvious  consideration  can  bHnadr  at  each  fixed 
w^elength,  the  pdariiation  signal  is  smaDer  in  photosphere,  and  stronger  in  the  umbra,  and  t^ 
m.^cls  are  very  different  each  other;  the  signal  inaeases  with  the  mcKnation,  both  for  longitudinal 
and  transverse  field;  there  is  quite  a  big  difference  between  wavelengths  distad't  as  low  as^20  mA 

«e  of^% ^  linearity 

are  of  2-3%  max,  the  s^e  is  not  completely  true  for  By,  where  these  deviations  can  reach  10-12% 

If  we  assume  a  power  law  relationship  between  Bj  and  P<j: 


Bt=C3P‘* 


(8) 


the  best  exponent  for  Pq  would  be  0.55-0.58,  instead  of  0.5. 

Brides  the  validity  of  the  WFA  errors  typical  of  this  calibration  method  are  related  to  velocity 
effects  line-of-sight  inclination  dependence,  model  dependence.  Trying  to  quantify  these  errors 
we  will  assuiM  as  rnean  cabbration  curves  for  our  magnetograms  the  ones  corresp^ding  to  the 
penumbral  model,  at  A=6122.04  A  and  7  =  30*  (for  Bi;  y  =  60“  for  Bj>).  ^ 

mV ‘-nefi-nafron  dependence.  Our  calibration  curve  corresponds  to  a  particular  value 
f  7,  but  nobody  can  tell  us  a  prtort  what  the  real  inclination  of  the  field  U.  Animate  of  the 
error  related  to  thu  assumption  is  given  in  Figs.  20a  and  21a,  for  Bl  and  Bj  respectively  The 
'"'"'T  polarisation  signals  corresponding  to  different  7  valu«,  for  the  same 

wavelength  and  model  and  converted  into  magnetic  fields  using  the  mean  curve  described  above 

Fe/ority  effecU.  A  velocity  field  will  shift  the  spectral  fine  of  a  quantity  AA  =  vX/c  and 
therefore  reduce  the  valiity  of  the  mean  calibration  curve,  since  we  will  be  at  a  dLent  wavelenSh 
from  the  working  one  Figs.  20b  and  21b  show  the  effects  of  a  velocity  field  of  ±lKm/!«  ^fa 
tjTical  value  for  solar  f^turw)  that  translates  mto  a  shift  of  ±20  mA.  As  in  the  previous  «ample 

the  curro  are  computed  usmg  the  polarisation  corresponding  to  the  adjacent  wavelenrths  fnd 
converted  into  magnetic  fields  through  the  mean  curve  wavelengths,  and 

The  errors  are  similar  b  both  the  and  By  cases,  and  of  the  order  of  15%. 

J.d.5  Model  dependence.  Usually,  many  different  features  are  present  b  the  FOV  therefore  the 

the  polarisation  signal  conespondmg  to  A=6122.04  A  and  7  =  30*  f60*t  for  the  t>,r J  j  1 

and  converted  it  through  the  mean  curves.  The  result.  L  sho^  b  Pb  2^  ’ 

inhe'fHd'^  “  *  Pbotospheric  situation  btroduci  an  undLttelb 

field  the  umbral  case  we  have  an  overestimate  of -15-20%.  In  the  transverse 

field  we  have  a  similar  overestimate,  of  ~  15-20%  b  both  the  situations. 

field^tbTd  ’  described  will  compensate  each  other,  like  the  case  of  umbral 

field  with  downward  motion,  but  of  course  we  cannot  count  on  this  kbd  of  "tricks-,  especiaUy 

*Note  tl3*t  below  -V  200  G  the  detemiinition  of  the  field  u  very 
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b€cau8c  we  don’t  know  anything  about  the  velocity  or  other  thermodynamical  parameter!  of  the 
•olar  feature!.  TVying  to  be  peMimietic,  we  built  the  i<M:alled  •everything  wrong*  curve,  where 
we  take  the  linear  polariiation  for  the  umbral  model,  at  >=8122.02  A  and  7  =  00*,  and  calibrate 
it  with  the  mean  curve:  the  rctiilt  u  in  Pig.  23,  and  the  enon  can  be  m  big  aa  50%. 

4.  DERIVATIVE  METHOD  IN  WEAK  FIELD  APPROXIMATION. 


4.a  The  ’’theory". 

In  a  recent  paper,  Jefferie!  and  Mickey  (1991)  extended  the  range  of  applicabffity  of  the  weak 
field  approximation,  aa  defined  in  Sect.  S.a,  provided  that  the  derivative!  <rf  the  Voigt  function 
describing  the  line  profile  fall  off  fast  enough.  At  the  very  end,  thii  tranalatea  into  the  posaibility 
of  using  the  WPA  for  value!  of  »/fl  aa  big  aa  1.5,  maintaining  a  predaion  of  15-20%,  if  we  work  at 
about  -  3  Ain  from  the  line  center;  this  result  ia  rather  mdependent  from  the  adopted  atmospheric 
model. 

The  equations  that  describe  this  "new*  WFA  are  (JM): 


Rcoey  = 


>1  V 
91  dl/dX 


(7) 


B8in7 


_  k»  /A>(g»  -I-  ir»)VJ 

9lS 


(8) 


where  dl/dX  ia  the  alope  of  the  intensity  profile,  A>  is  the  of&et  from  line  center.  ki  =  5.4  x  10*  x 
(6300/>)’,  4j  =  8.2  X  10^  x  (6300/>)*,  with  wavelengths  measured  in  A  and  B  in  Gauss. 

Eqa.  (7)  and  (8)  are  valid  for  all  the  linea,  and,  at  least  in  principle,  the  method  shouldn’t  be 
affected  by  model  dependence  or  velocity  effects  (as  far  as  the  shift  is  small  enough  to  keep  us  in 
the  range  of  validity  of  the  WFA).  This  seems  therefore  a  very  attractive  method  for  the  calibration 
of  imaging  magnetographs:  the  only  requirement  is  the  acquisition  of  2  or  more  images  close  in 
wavelength,  so  to  estimate  the  value  otdI/dX. 


An  obvious  first  obstacle  is  that  at  3  (or  more)  AA^,  from  the  line  center,  the  polarisation 
signal  18  extremely  small,  and  hence  very  difficult  to  detect:  b  Fig.  19,  e.g.,  we  see  that  at  -2£^Xo 
(>=6121.96  A)  the  linear  polarisation  reaches  detectable  values  only  for  B  >  500  G-  at  3A>n  this 
lower  limit  bereases  to  ~  800  G.  ’ 

We  c^  check  if  a  sort  of -compromise*  between  the  distance  from  Ibe  center  and  value  of  fractional 
polarisation  will  work,  may  be  at  the  prices  of  a  bigger  uncertabty  on  the  field  values. 

Fig.  24  shows  the  curves  obtained  with  Eq.  (7)  for  the  same  model,  field  values,  bclinations  and 
wavelengths  of  Fig.  8.  The  slope  of  the  mtensity  profile  has  been  evaluated  with  the  ■analytical- 
derivative  of  the  Stokes  I;  the  spectral  smearing  due  to  FPI+blocker  hasn’t  been  taken  bto  account 
Fig.  25  shows  the  correspondent  percentage  errors,  defined  as  (Bi„,  -  B„,)  x  100/J5i„  .  The 
agreement  between  the  bput  and  the  computed  fields  is  much  better  further  b  the  wing  (as 
expected)  and  for  smaller  fields,  and  differences  due  to  the  bclbation  values  still  exist.  In  the 
range  of  wavelengths  previously  defined  on  the  base  of  the  polarisation  signals,  the  errors  that  we 
get  using  eq.  (7)  are  withb  10%  up  to  a  1000  G,  but  get  much  worse  for  bigger  fields. 

Figs.  26-27  Me  the  same,  for  the  transverse  field  case,  i.e.,  nsbg  eq.  (8).  The  genwal  trend  is 
the  same  as  before,  but  the  amplitude  of  the  errors  is  bigger.  For  the  case  of  pure  transverse  field 
(7  -  90  )  the  errors  are  quite  small  everywhere,  but  this  is  unfortunately  not  true  b  the  other  2 
cases,  even  though  we  remab  b  the  range  of  ~20%  foreseen  by  JM. 

In  the  ;  otospheric  case,  the  situation  u  better  for  the  transverse  field,  and  a  little  worse  for  the 

longitudmal;  the  contrary  u  true  m  the  umbral  model,  but,  agab,  we  remab  b  the  20%  of  error 
regime. 

Note  that  the  errors  are  all  positive,  i.e.,  the  derivative  caUbration  gives  a  systematic  underestimate 
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4.b  EfTects  of  spectral  smearing.  FPl 

No  mention  of  the  effects  of  spectral  smearing  is  given  in  JM  but,  as  we  saw  before,  the  presence 
of  a  smearing  function  introduces  significant  changes  respect  to  the  original  situation. 

Figs.  28  and  29  are  the  correspondent  of  25  and  27,  for  profiles  convolved  with  the  FPI  transmission 
profile  {FWHM=150  mA). 

For  Bi  we  see  that  there  is  a  smaller  error  for  the  case  of  pure  lonptudinal  field  at  all  wavelengths, 
aad  that  the  convolution  has  somehow  reduced  the  "wild"  behaviors  m  the  wavelengths  close  to 
line  center,  but  the  differences  among  inclinations  are  bigger  than  before;  we  saw  in  Sect.  3.b  that 
the  ratio  V/I  is  more  homogeneous  for  the  convolved  profiles  respect  to  the  original  ones,  so  the 
cause  of  this  discrepancy  must  belong  to  the  factor  dl/dX. 

Indeed,  the  derivative  of  the  convolved  intensity  profiles  is  strongly  dependent  on  the  7  values; 
something  similar  to  the  effect  described  in  Sect.  S.b  about  the  fractional  linear  polarisation  is 
happening  here  too:  the  convolution  process  "brings"  part  of  the  signal  at  line  center,  extremely 
dependent  on  the  field  and  the  inclination,  to  the  wings,  causing  the  behavior  shown  in  Fig.  28. 
So,  the  l^t  that  we  can  say  is  that,  in  the  range  of  wavelengths  6122.0-6122.08  A,  the  errors  for 
Bl  rem^  within  the  20%. 

The  transverse  field  (Fig.  29)  shows  tremendously  enhanced  errors  and  a  very  large  spread  within 
different  inclinations. 

The  reason  for  the  increase  m  the  error  curves  is  easily  understood  looking  at  formula  (8):  the 
proportionality  between  the  linear  polarisation  and  the  slope  of  the  intensity  profile  is  ruled  not  by 
the  Bt  only,  but  also  by  the  numerical  factor  AX,  that,  being  the  offset  from  line  center,  doesn’t 
change  either  we  convolve  the  profiles  or  not.  In  the  un-smeared  profiles,  the  shapes  of  Q  (U) 
and  dl/dX  are  steep  enough  to  compensate  for  the  the  AX  factor,  while  this  is  not  true  for  the 
convolved  profiles.  Fig.  30  visualises  the  point:  using  the  original  U,  Q,  I,  the  quantity  defined 
by  (8)  asymptotically  tends  to  the  input  field;  using  the  convolved  ones,  the  trend  that  we  get  is 
basically  the  shape  of  AA. 

We  need,  therefore,  an  "effective"  AA,  somehow  weighted  with  the  spectral  transmission  profile  of 
the  FPI®,  but  of  course  not  dependent  on  the  field  strength  or  the  inclination  value. 

The  best  way  we  could  think  to  determine  this  AA,/y  was  just  a  triaj-and-error  method,  trying  to 
minimiie  the  enor  curves  of  Fig.  29.;  for  the  pcnumbral  and  umbra]  models,  ancf  wavelengths  in 
the  interval  6122.0-6122.06  A,  the  best  value  of  AA.yy  u  ~  0.65  x  AA,  while  for  the  photospheric 
model  is  0.5  X  AA.  *  ^ 

Fig.  31  shows  these  "new”  error  curves,  that  now  start  and  end  within  ±30%;  however,  the  big 
dependence  on  the  7  values  is  still  present,  a  little  enhanced  respect  to  the  un-convolved  curves. 
4.C  Effects  of  spectral  smearing.  Blocker 

As  described  in  Sect.  3.c,  the  net  effect  of  the  blocker  u  to  reduce  the  intensity  transmission 
of  an  appreciable  quantity,  even  if  positioned  at  the  working  wavelength,  whUe  the  polarization 
signals  are  less  affected. 

This  implies  that  the  factor  dl/dX  will  be  smaUcr  and,  using  the  derivative  calibration,  the  output 
field  will  be  an  overestimate  of  the  true  one.  The  enor  curves  of  Fig.  28  will  therefore  shift  towards 
more  negative  values,  even  though  they  will  keep  the  same  shape;  the  amount  of  the  shift  can  be 
as  large  as  10%. 

The  shift  for  Bj  is  limited  to  ~3%,  since  the  factor  dl/dX  is  under  square  root  (Eq.  (8)).  Fig.  32 
show  the  same  enor  curves  of  Figs.  28  and  31,  for  B^,  and  IJ,-  and  a  smafler  range  of  wavelengths 
considering  the  blocker  factor.  ’ 

Another  possible  source  of  enors  is  the  misplacing  of  the  blocker.  The  transmission  of  the 
blocker  falls  off  quite  rapidly,  respect  to  the  width  of  the  fine,  and  therefore  a  misplacing  of  even 
few  tenths  of  mA  respect  to  the  working  wavelength  causes  noticeable  differences.  The  net  result 

conToIution  won’t  be  u»eful,  .ince  the  operator  convolution  doeanH  commute  with  the  operator 
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depends  on  the  nmount  end  iirteUon  of  the  blocker  shift:  Pig.  S3  shows  fm  the  same  format  as 
Pig.  32)  the  error  ciurves  for  a  shift  of  the  blocker  of  40  mA,  Pig.  34  for  a  shift  of  ^  mX. 

A  precise  tuning  of  the  blocker  is  hence  an  absolute  necessity  for  this  calibration  technique. 

4.d  Effects  of  finite  difference 

All  the  previous  curves  have  been  computed  using  the  "analytic*  derivative  of  the  profiles; 
with  the  data  we  have  to  use  a  finite  difference  between  intensities  at  2  different  wavelengths 
positions  in  the  wing  of  the  line.  Here  we  will  consider  the  effects  of  doing  so,  but  for  profiles 
convolved  with  the  FPI  alone. 

The  estimate  of  the  derivative  will  be  close  to  the  real  value  if  the  Kparation  in  wavelength  of  the 
intensitiM  images  is  small,  but  on  the  other  hand  if  the  images  are  too  close,  the  corresponding 
signal  will  be  too  low  to  detect.  Considering  the  convolved  intensity  profiles,  an  interval  of  20  mA 
will  give  a  d//f  ~  2-3%,  that  we  can  detect. 

Given  the  working  wavelength  A<,  we  can  estimate  the  derivative  either  as  A/  =  /^^  -  /j.  ©r 
A/  =  -  fiii  In  the  first  case  we  have  a  value  bigger  than  the  real  one,  m  the  second  a’^ue 

smaller.  This  translates  into  a  shift  of  the  error  curves  <5%,  as  shown  in  Kg.  35.  As  explained  in 
Sect.  4.C,  the  presence  of  the  blocker  mtroduces  an  overestimate  of  the  field  that  can  be  as  big  as 
10%  m  the  Bi  case;  to  avoid  that  the  error  introduced  by  the  finite  difference  sums  with  this  we 
need  to  use  the  A/  =  as  estimate  of  the  derivative. 

5.  DETERMINATION  OF  THE  AZIMUTH 
The  asimuth  of  the  field,  x,  «  usually  determined  through  the  formula: 

X  =  0.5  arctan(£7/Q) 

Its  value  should  be  wavelength-independent,  as  far  as  we  remain  in  the  line  wings.  Fig  36  shows 
the  values  of  the  asimuth  (input  x  =  0*)  determined  with  eq.  (9),  for  the  same  range  of  wavelengths 
of  the  previous  figures.  For  the  case  of  pure  transverse  field  (7  =  0*)  the  output  value  is  perfectly 
sero,  but  for  the  other  inclinations,  and  for  small  fields,  there  is  an  error  <  10*.  At  bigger  fields 
the  error  tends  asympotically  to  sero;  for  wavelengths  closer  to  the  line  center,  though,  the  errors 
increase  ag^.  ’ 

When  we  convolve  with  the  usual  Fabry-Perot  profile  of  150  mA  FWHM,  the  curves  become  more 
homogeneous,  and  we  have  an  undetermination  of  10  -  12*,  for  aU  the  wavelengths,  up  to  800  G- 
over  this  value  the  errors  reduce  to  ~  5®  (Fig.  37).  ’ 

6.  SUMMARY  AND  CONCLUSIONS 

As  mentioned  in  the  Introduction,  the  conversion  of  polarisation  into  magnetic  fields  is  a  real 
difficult  task  and,  indeed,  what  we  saw  in  the  previous  sections  is  not  very  encouraging. 

We  anJys^  two  methods  that  could  be  suitable  for  the  calibration  of  Vector  Magnetographs, 
the  derivative  method  and  the  calibration  curve  method,  and  commented  the  determinationof  the 
asimuth  of  the  field.  Here  we  will  try  to  summarise  the  resulU,  and  pve  some  suggestions  for 
further  improvements. 


Se  Atene  pianye  ... 

Derivative  method  The  derivative  calibration  method  works  reasonably  well  if  we  go  in  the  far 
wings  of  the  Une;  here,  the  average  errors  Ue  within  15-20%,  both  for  longitudinal  and  transverse 
fields,  up  to  1500-2000  G  (Figs.  25  and  27).  However,  this  is  not  a  general  stetement,  since  the 
range  of  validity  greatly  changes  from  line  to  line. 

When  we  introduce  a  convolution  factor,  as  in  the  case  of  a  filter  magnetograph,  things  change 
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drastiodly:  the  error*  «e  •omehow  itUI  accepteble  m  the  Bi  caw,  but  for  By  we  go  «P  to  30% 
or  more,  and  there  is  a  strong  dependence  of  thew  error*  to  the  bclinatlon  of  the  field,  7.  A  filter 
with  a  narrower  passband  respect  the  one  actually  present  b  the  SVMG  would  reduce  this  effects. 
The  presence  of  the  blocker  is  another  source  of  error  to  consider  using  this  method  (Fig*. 
i4).  The  two  major  problem*  are  its  very  narrow  transmission  profile  respect  to  the  width  of  the 
line  and  it*  not  completely  reliable  tunabiUty;  a  tunahle  blocker  with  a  bigger  passband  1.6 

A  FWHM,  but  not  more,  sbce  the  a4jacent  order*  of  the  FPI  lie  at  -  3.3  A  apart)  would  much 
improve  the  situation. 

-ne  method  U  not  much  affected  by  velocity  effect*  or  by  -model  dependence-,  b  the  senw  that 
the  errors  are  comparable  in  the  range  of  wavelength  6122.0-fll22.0«  A  and  for  the  three  models 
The  neewity  of  using  finite  differences  instead  of -analytical-  derivatives  doesn’t  ’introduce  error* 
too  big  (<  5%).  if  we  limit  the  step  to  20-30  mA,  since  the  Ibe  is  wide  enough  to  change  very 
smoothly  in  this  interval.  The  same  wouldnH  be  true  with  a  narrower  fine,  like  the  Fel  5250  (see 
Caussi  ei  al.,  1991).  Sbce  we  are  subtracting  image*  acquired  at  different  wavelengths,  though, 
we  must  pay  special  attention  to  some  factors,  like  eventual  dependence  on  wavelength  of  the 
transmission  of  the  instrument,  differences  b  the  transparency  of  the  sky,  etc.  The  derivative 
factor  seems  to  be  extremely  sensitive  to  misalignmcnU  of  the  images,  and  this  can  be  a  wriou* 
problem  m  the  actual  version  of  the  SVMG,  sbce  the  btensity  images  are  constructed  on  a  quite 
long  time  lap  (between  20  sec,  b  the  -old-  mode,  and  20  minutes,  b  the  -new*  one),  with  the 
unavoidable  motion  of  the  features  in  the  FOV.  A  smooth  of  the  btensity  image*  before  the 
subtraction  will  improve  the  situation. 

b  summary,  b  the  actual  version  of  the  SVMG,  the  derivative  method  can  give  u*  results  with  an 
mdetermination  of  30%  or  worse. 


...  Sparta  non  ride 


Calibration  curve  methol  As  seen  b  Sect.  3.a,  the  relationships  between  fractional  polarisation 
and  magnetic  fields  are  not  completely  Ibear,  even  after  the  convolution  with  a  filter  150  mA  wide 
and  therefore  the  use  of  Eqs.  (3)  and  (4)  as  calibration  curves  gives  us  enors  that  can  be  up  to 

12%  for  the  ^  case.  However,  we  could  solve  the  issue  using  different  curves,  like  a  power  law  for 
the  linear  polansation. 

•^e  convolution  with  the  FPI  profile  increases  the  linearity  of  the  curves,  but  has  ^o  opposite 
effects  on  the  differences  due  to  the  inclination  values  for  Bj,  and  Bx  (Figs.  11  ^d  12).  A 
fflbr  with  a  narrower  passband  would  be  better  for  the  transverse  field,  even  though  quite  big 
<^erenccs  we  present  b  the  original  profiles,  but  would  reduce  the  homogeneity  of  the  curves  for 
the  longitudinal  field,  and  vice  versa. 

The  blocker  doesn’t  affect  much  the  precision  of  the  method. 


Errors  intrinsic  b  the  method,  as  described  b  sect.  3.d,  are  due  to  velocity  and  bcibation  effect* 
and  Ic^  to  a  total  enor  that  could  be  as  big  as  30%;  an  accurate  choice  of  the  calibration  curve 
could  however  rmnimise  the  effects  due  to  the  7  value*. 

The  problems  btroduced  by  the  choice  of  a  particular  model  could  be  overcome  adopting  the  right 
model  for  each  feature  in  the  FOV,  even  if  this  wiH  heavily  affect  the  computational  time;  b  L 
case  the  ^certainties  of  the  models  themselves,  or  b  the  determbation  of  which  model  to  use  for 
a  particular  pomt,  will  remab. 


b  sui^ary,  once  we  t^t  and  adopt  a  particular  atmospheric  model,  the  calibration  curve  method 
can  give  us  results  with  an  bdetermbation  (btemal  consistence?)  of  —  20%  This  result  u 
not  strongly  affected  by  the  characteristic*  of  the  SVMG,  but  will  equaUy  apply  to  any  filter 


Both  the  methods  will  suffer  by  the  presence  of  scattered  light;  as  an  example,  the  MarshaU 
group  correc^  for  scattered  light  comparing  their  results  with  the  values  of  Mt.  Wilson  (obtabed 
spectrographically),  and  the  correction  factor  is  around  8  (see  Hagyard  et  al.,  1988,  App.  E). 
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All  thU  coDiideratioiu,  plu*  the  comment*  c(  Sect.  6,  lead  ui  to  eome  conehuion*: 

1.  The  derivative  method  wiU  work  better  with  a  very  narrow  filter  (we  think  that  the  method 
would  be  very  luiUble  for  «pectrographic  obeervationi)  and  need*  a  tunable  blocker  to  that 
the  global  trarumisaion  profile  U  the  one  of  the  filter  only.  An  additional  improvement  would 
come  by  apeeding  up  the  acquiaition  rate.  The  method  u,  however,  not  suitable  for  all  the 
linear  amall  AAp  value*  will  greatly  reduce  the  range  of  vJidity  of  the  WPA. 

2.  If  we  choose  to  use  the  calibration  curve  method,  the  filter  shouldnH  have  a  paaaband  smaUer 
than  ~  100  mA,  or  there  will  be  too  much  room  for  errors  related  to  velocity  fields  (Hagyard 
et  al.,  1988,  Chap.  5);  over  this  value,  there  are  equivalent  pro*  and  cons,  as  reminded 
above.  The  method  is  most  suitable  for  line*  with  a  high  gi  factor,  since  thi*  wiU  improve 
the  polarisation  signal. 

3.  In  any  case,  the  vector  magnetic  fields  that  we  can  determine  with  filter  magnetograph*  are, 
at  the  best,  a  very  rough  estimate  of  the  real  quantities;  we  can  expect  errors  of  2-500  G  in 
the  transverse  field,  up  to  200  G  for  the  longitudinal,  and  ~  10*  for  the  asimuth.  This  is 
based  only  on  the  analysis  of  the  calibration  methods,  i.e.  without  considering  other  sources 
of  noise,  Uke  instrumental  polarisation,  scattered  light,  mhomogeneities  in  the  transmission, 
etc.  Nevertheleu,  a  filter  magnetograph  can  give  us  information  about  the  whole  structure’ 
of  an  active  region,  with  a  good  spatial  resolution  and,  in  general,  we’ll  be  able  to  follow  the 
temporal  evolution  of  the  structures  in  considerable  detail. 
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FIGURE  CAPTIONS 


Fig.  la.  SVMG  Etalon  ET75  transmission  profile.  The  transmission  has  been  arbitrarily 
•et  to  50%.  ^ 

interference  blocker,  at  T-50*C  and  titt=0*.  The  transmission 

value  18  ~27%. 

^  (solid)  and  V  (dashed)  for  Cal  6122,  photospheric  model 
and  3  different  sets  of  magnetic  parameters.  The  max.  circular  polarization  is  --20%.  ’ 

Fig.  S.  As  Fig.S,  but  for  Stokes  Q  (soHd)  and  U  (dashed).  The  maximum  signal  is  10% 
for  the  case  of  7  =  60*.  ’ 

Figs.  4~5-  As  Figs.  2-S,  for  Moe-Maltby  penumbral  model. 

Figs.  6-7.  As  Figs.  S-S,  for  Maltby  umbral  model. 

8.  Fractional  circular  polarization  vs  input  longitudinal  B,  for  the  penumbral  model 

The  signal  increases  towards  Une  center,  but  the  same  do  magnetooptical  effects  ^  de¬ 
pendence  etc.  '  ' 

9.  Fractional  linear  polarization  vs  input  transverse  B,  for  the  penumbral  model 
^e  same  considerations  of  Fig.  8  apply.  The  signal  is  extremely  smaU  for  fields  <  500 


Fig.  10.  The  s^e  of  Fig.  P,  but  on  logarithmic  scale.  A  value  of  linear  polarization  of 
10  *  corresponds  to  1-200  G. 

F^  11.  Effect  of  the  FPI  (spectral  selector)  on  the  Stokes  profiles.  The  curve  are 
shaUower  than  the  original  ones,  but  the  values  of  fractional  polarization  are  biwer  in 

the  far  wings  of  the  line.  The  details  of  Q  and  U  at  line  center  have  been  wasW  out 
almost  completely. 

smearing  intro- 

duced  by  the  Fabry-Perot.  There  is  a  strong  7  dependence  in  the  linear  polariaation. 

fX  M  ‘’’'i'"”-  T"‘  inclinations 

(7  -  90  and  30  ).  The  signal  is  similar  m  the  far  wings,  b)  The  same  as  before  but 

involved  with  the  FPI  transmission  profile  (FWHM=150  mA).  The  signal  even  b  the 
tar  wings  is  very  different. 

Fig.  16.  Effects  of  the  blocker  prefilter  on  the  Stokes  parameters. 

of  I2-IS-I4,  but  introducing  the  blocker  factor.  These  are  the 
curves  tiid;  v?ill  be  used  as  ccuih7Xt,tioTi  curves. 

Fig.  eo.  a)  B  long,  obtained  from  drcular  polarization  correspondent  to  3  different 
lme-of-«ght  incbnations  and  calibrated  with  the  mean  curve,  »a  the  true  B  long.  The 
^rs  he  within  10%.  b)  B  ong.  obtmned  from  circular  polarization  conespondent  to  3 

caUbrated  with  the  mean  curve,  vs.  the  true  B  long.  The  errors 


F^.  21.  a)  B  tran.  obtained  from  linear  polarization  correspondent  to  3  different  line- 
of-8ight  inclinations  and  calibrated  with  the  mean  curve,  vs,  the  true  B  tran  The  errors 
are  pretty  big,  with  a  maximum  of  ~  40%  for  7  =  90“.  b)  B  tran.  obtained  from  linear 
polanzation  correspondent  to  3  different  wavelengths  and  calibrated  with  the  mean  curve 
VB.  the  true  B  tran.  The  errors  are  similar  to  those  of  Fig.  20  b).  * 


ru  *  wrong  model.  Umbral,  penumbral  and  photospheric  signal 

^brated  with  a  penumbral  curve,  for  the  case  of  B  long.  The  errors  can  go  up  to  a  30% 

"’Snal-  b)  Same  as  a),  for  the  transverse  field.  The  errors  remain 

within  16-20%. 


Fig  2S.  Mi^ng  of  different  factors:  wrong  model,  redshifl  of  1  Km/sec,  different  inch- 
nation,  for  the  case.  The  total  error  is  around  50%. 

Fig.  24.  B  longitudinal  computed  with  the  derivative  method  (for  penumbral  model  3  -y 
and  8  wavelengths  as  before)  vs.  input  B  long.  No  spectral  smearing. 

Fig.  25.  Percentage  error  curves  correspondent  to  the  curves  of  Fig.  24.  Far  in  the  wines 
the  errors  are  within  ~  10%,  ®  ’ 


Figs.  26-27.  Analogous  of  Figs.  24-25,  for  the  transverse  field.  The  errors  are  worse 
this  case,  especially  for  7  =  30  and  60®. 


in 


Fig.  28.  As  Fig.  25,  considering  the  spectral  smearing  introduced  by  the  FPI  There  is  a 
bigger  spread  among  inclinations  that  before. 

29.  Like  Fig.  S8,  hnt  for  Bj.  The  errors  are  tremendously  increased,  due  to  the 

m  forrrmla  fSl  in  ’ 


Fig. 


factor  AA  in  formula  (8)  in  the  text. 

Fig.  SO.  Spectral  behavior  of  the  Bt  computed  with  the  derivative  method  for  un¬ 
smeared  (solid)  and  convolved  (dashed)  profiles.  In  the  first  case  the  output  tends  to  the 
true  value,  in  the  second  the  shape  of  AA  dominates. 

=  O  MAA.  Tht  erron  remain  within 
dzoU/o,  but  the  spread  between  inclinations  is  very  big. 

F^.  S2.  MuenM  of  the  blocker  in  the  derivative  method.  The  global  effect  is  to  shift 
the  curves  towards  more  negative  values  (of  about  10%  for  B^  and  3%  for  Bt). 

Fics.  SS-Sj.  Effects  of  blocker  misplacing.  In  the  first  figure,  the  blocker  is  set  at  -f-40 
mA  respect  to  the  working  wavelength;  in  the  second,  to  -40  inA. 

Fig.  55.  Effects  of  fmte  difference.  The  estimate  of  dlfdX  with  a  finite  difference  between 

StlTht  m  od y  H* 

determination  for  un-convolved  profiles.  The  error  remains  <  10®  in 
all  the  situations.  _  « 


Fig.  57.  Same  as  36,  but  for  convolved  profiles  (FPI  only).  The  errors  are  a  little  increased 
lor  smaller  fields,  but  ~  10®  again 
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Table  1.  ci  values  for  3  models,  3  inclinations  and  4  wavelengths 


XX{Angstrom3) 


6122.00 

6122.02 

6122.04 

6122.06 


Photosphere 
7  =  0,30,60* 

24520,24315,23294 
20197, 19996, 19068 
16650,16450,15615 
13840, 13660, 12920 


Penumbra 
7  =  0,30,60* 

16390,16540,15688 

14255,14410,13600 

12490,12668,11900 

11110,11335,10630 


Umbra 
7  =  0,30,60* 

12720,12670,12240 
11840,11780,11320 
10950, 10880, 10380 
10080,10000,9470 


Table  2.  C2  values  for  3  models,  3  inclinations  and  4  wavelengths 


XX{Angstroms) 
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12300,13080,10530 

11660,11690,9420 

10530,10540,8500 

9610,9650,7790 


Umbra 
7  =  30,60,90* 

11860,12115,10230 

10800,11040,9385 

9800,10020,8580 

8870,9081,7846 
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